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Abstract We study gravitational radiation reaction in the equations of motion for 
binary systems with spin-orbit coupling, at order (i>/c) 7 beyond Newtonian gravity, 
or 0(v/c) 2 beyond the leading radiation reaction effects for non-spinning bodies. We 
use expressions for the energy and angular momentum flux at infinity that include 
spin-orbit corrections, together with an assumption of energy and angular momentum 
balance, to derive equations of motion that are valid for general orbits and for a class 
of coordinate gauges. We show that the equations of motion are compatible with those 
derived earlier by a direct calculation. 
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1 Introduction and Summary 

The backreaction of the emission of gravitational radiation on the system emitting 
the radiation is a problem both of formal interest within general relativity and of 
practical interest for gravitational-wave detection. A leading candidate source for laser- 
interferometric gravitational-wave observatories, both on the ground and in space, is the 
radiation-reaction induced inspiral of a binary system of two compact objects (black 
holes or neutron stars). In order to develop accurate theoretical predictions for the 
gravitational waveforms emitted by such systems, one must know their evolution under 
the dissipative effects of gravitational- wave emission to high accuracy. In addition, 
particularly for systems containing black holes, the effects of spin may be important. 
Spin-orbit and spin-spin couplings can result in precessions of the bodies' spins and of 
the orbital angular momentum, leading to modulations in the gravitational waveform 
P], and can affect the rate of decay of the orbit 2,3 . 
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As a result, substantial effort has gone into determining the effects of spin in binary 
systems. Except for the final few orbits, much of the inspiral of such systems can be 
described by the post-Newtonian approximation, which is an expansion of Einstein's 
equations in powers of e ~ (v/c) 2 ~ Gm/rc 2 , where v, m and r represent typical 
velocities, masses and separations in the system, and G and c are the gravitational 
constant and speed of light. Each power of e represents one "post-Newtonian" (PN) 
order in the series (e 1 / 2 represents one-half, or 0.5PN orders). 

Formally, spin effects first enter the equations of motion at the 1PN level, and 
have been derived by numerous authors from a variety of points of view, ranging from 
formal developments of the GR equations of motion in multipole expansions [4j[5] , to 
post-Newtonian calculations [6], to treatments of linearized GR as a spin- two quantum 
theory 7,8 . For a review of these various approaches, see [5]. 

Spin also affects gravitational radiation reaction, and radiation reaction can affect 
spin; it is straightforward to show that such effects first occur at 3.5PN order. In 
earlier work, we derived, from first principles, the radiation-reaction effects of spin-orbit 
and spin-spin coupling, by integrating the post-Newtonian hydrodynamic equations of 
motion, including 1PN, 2.5PN and 3.5PN terms, over bodies consisting of rotating fluid 
10, 11]. As a check, we found that the loss of energy and angular momentum (including 
spin) induced by the radiation reaction terms, matched precisely the expressions for 
energy and angular momemtum flux derived by Kidder et aZ.(2][3]- 

An alternative approach to obtaining equations of motion with radiation reaction 
at higher PN orders was studied by Iyer and Will jl 2 1 1 13 j - There, we wrote down the 
most general form that the 2.5PN and 3.5PN radiation-reaction terms could take in 
the equations of motion for a binary system of spinless bodies, in terms of arbitrary 
coefficients. We then used the assumption of energy and angular momentum balance, 
combined with energy and angular-momentum flux expressions accurate to PN order 
beyond the quadrupole approximation to impose constraints on the arbitrary coef- 
ficients used in the equations of motion. After taking into account a fundamental 
ambiguity in the definitions of energy and angular momentum at 2.5PN and 3.5PN 
orders, we were left with equations of motion with coefficients that are fixed up to 
two arbitrary coefficients at 2.5PN order and 6 arbitrary coefficients at 3.5PN order. 
It was then straightforward to show that these eight degrees of freedom correspond 
precisely to the effects, mapped onto the two-body equations of motion, of coordinate 
transformations at the relevant PN orders. At 2.5PN order, for example, one choice of 
the two arbitrary coefficients gives the equations in the so-called Burke-Thorne gauge 
[14] . in which the radiation reaction terms are obtained from a gradient of the potential 
(G /5c 2 )x t d 5 / <y> /dt' , where J <4:,> is the trace-free moment of inertia tensor of the 
system, while another choice gives the D amour- Deruelle gauge, which is more directly 
tied to harmonic gauge [151116] . For spinless systems, this approach was extended to 
determine the 4.5PN terms in the equations of motion using flux expressions accurate 
to 2PN order beyond quadrupole [17] . 

It is the purpose of this paper to extend this approach to include spin-orbit radiation 
reaction effects at 3.5PN order. We assume that the equation of motion for the relative 
vector x = xi — X2 in a binary system may be written in a PN expansion in the form 

m . . 

a = --jn + a PN _ so + • • . + a 2 . 5 PN + • • • + a 3 . 5 p N _ S o , (1) 

where r = |x|, n = x/r and m = mi + 7T12; apN_go is the post-Newtonian spin-orbit 
contribution, a2.5PN is the leading radiation-reaction contribution, and a3.5PN_so is 
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the 3.5PN spin-orbit contribution. Here and for the rest of this paper, we use units in 
which G = c = 1. We have not displayed the point mass 1PN, 2PN, 3PN, and 3.5PN 
terms, as they will play no role in our analysis. We also will not use the bookeeping 
parameter e explicitly to keep track of PN orders, since we will be considering only 
specific orders. It is sufficient to recall that, since spin scales as rnvr, then S/r 2 ~ 
v(m/r) ~ e 3 / 2 . This, plus explicit labelling of terms throughout, should make clear the 
PN order of terms being discussed. 

We then write down the most general 3.5PN spin-orbit expression that (a) con- 
tains terms each involving a single spin (either Si or S2), (b) is a vector, and (c) is 
antisymmetric under the interchange 1^2. This turns out to involve 30 arbitrary 
coefficients. 

Because the bodies have intrinsic spin, we must make an assumption about their 
spin evolution. At 1PN order, we assume that they obey the standard spin-orbit pre- 
cession equations (see Eq. ([2)| below). These 1PN equations produce only precession; 
the magnitudes of the spins do not change. At 3.5PN order, we likewise assume that 
gravitational radiation reaction produces only precessions of the spin. This is a rea- 
sonable assumption, because, for a rotating axisymmetric body, it is impossible to see 
how gravitational radiation can cause it to spin up or down, to the 3.5PN order being 
considered. Such effects must involve specific couplings of radiation to deformations 
of the bodies, either due to rotational flattening or due to tidal couplings, which are 
beyond the scope of our assumption of almost point-like bodies with spin. In this case, 
we can then show that the most general 3.5PN expression for the evolution of each spin 
that (a) is a pseudovector; (b) depends only on the spin itself and on orbital variables; 
and (c) is orthogonal to the spin, can in fact be written as a total time derivative of 
spin and orbital variables, which can then be absorbed into a meaningless 3.5PN-order 
correction to the definition of the spin. 

Consequently, we can calculate the loss of energy and angular momentum using 
only the parametrized equations of motion and the 1PN spin precession equations. 
There is no contribution to the evolution of the spins at 3.5PN order. However, we 
must incorporate the freedom to add arbitrary terms of 3.5PN spin-orbit order into 
the definitions of total energy and total angular momentum, just as in the spinless case. 
There are 6 such terms in E and 26 in J. Thus there is a total of 62 arbitrary coefficients 
to be determined. We then equate the time derivative of these expressions for E and 
J with the corresponding expressions obtained from the far-zone gravitational-wave 
flux, including spin-orbit terms, as calculated by Kidder el al. an d compare them 
term by term. This leads to 54 constraints on the coefficients; however 4 of these 
constraints are not linearly independent of others, and thus we have 50 constraints 
on 62 coefficients, leaving 12 undetermined coefficients. Finally we show that these 12 
free coefficients in the equation of motion correspond precisely to the effects of 3.5PN 
order coordinate transformations, mapped onto the two-body equations of motion with 
spin-orbit coupling. 

The remainder of this paper provides details. In Sec. [2] we review the known equa- 
tions of motion and spin evolution through 2.5PN order. Section [3] applies energy and 
angular momentum balance to determine the 3.5PN spin-orbit terms in the two-body 
equations of motion, while Sec. [4] shows that the remaining undetermined coefficients 
are directly related to gauge freedom. Section [5] presents concluding remarks, while 
certain detailed formulae are relegated to Appendices. 
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2 Two-body equations of motion with spin-orbit coupling 

The PN-SO and 2.5PN terms in Eq. |TJ are given by conventional expressions 



apN-so = -U ™£n- (45 + 30- vx (45 + 3€) + 2fiix (4S + 3*) }, (2) 



&2.5PN = 



r 3 1 2 r 

8/j.m 
5r 3 



3(1 + (3)v 2 + I (23 + 6a - 9/3) — - 5/3r 2 
3 r 



(2 + a)u 2 + (2 - a)— - 3(1 + a)r 2 
r 



(3) 



where v = dx/dt is the relative velocity, /x = mim^/m is the reduced mass, 5 = >Si+52 
is the total spin, £ = {m2/mi)S\ + (mi/m,2)S2 is a second spin parameter, Ln = xxv 
is the orbital angular momentum per unit reduced mass, and r = v • n. 

The coefficients a and /3 in a 2 .5PN reflect the possibility of different gauges for 
expressing radiation reaction at 2.5PN order [121113] . The choice a = 4, f3 — 5 cor- 
responds to Burke-Thorne gauge [14], while the choice a = — 1, f3 = leads to the 
Damour-Deruelle radiation-reaction formula [151116] . Any choice of a and /3 leads to 
the same loss of energy and angular momentum at 2.5PN order, corresponding to 
quadrupole approximation energy and angular momentum flux. 

In defining spin, we must specify the center of mass of each body using a procedure 
commonly known as the "spin supplementary condition (SSC); the definition used in 
this paper corresponds to the value fcsSC = 1/2 (see [181131110] for further discussion). 

The equations of evolution for the spins may be written in the form 



5l = (5l) PN + ... + (5; 



U3.5PN-SO i 



(4) 



where 



(Sl)pN = ^-^N 



x Si 



3 m 2 \ 
2 mi) 



(5) 



and where the equations for spin 2 can be obtained by the interchange 1 ^ 2. We have 
not included conservative 2PN and 3PN contributions, and can show that the leading 
radiation reaction contributions come at 3.5PN order [10] . 

Up to 2.5PN order, the motion is conservative and the energy and angular momen- 
tum are constant. Including only the Newtonian and spin-orbit terms, they are given 
by 



E = E N = fj,\ -v 

2 r 



J = /xL N + S + — n x 

2r 



n x (45 + 3£) 



(6) 
(7) 



In our chosen spin supplementary condition, there is no spin-orbit contribution to 
the conserved energy, while J contains the orbital angular momentum, the total spin, 
and a PN spin-orbit contribution. These conserved quantities can be derived from the 
equations of motion by constructing ^dv 2 /dt = v • a, and d(x x v)/dt = x x a, and 
showing that, after substituting the equations of motion and spin- precession equations 
carried to the appropriate order, everything can be expressed as total time derivatives. 
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3 Spin-orbit radiation reaction via E and J balance 



We now write down the most general 3.5PN spin-orbit terms as 



a 3.5PN-SO 



5r 4 



As — (Ln ■ S) + B s -(l N ■ S) + C s fv xS + D s n x S 
r r 



+A £ — (L N ■$)+B^-(L N ■$) + C^r\ x$ + D^n x £ 



(8) 



The form of Eq. ([8| is dictated by the fact that it must be a correction to the Newtonian 
acceleration, (i.e. be proportional to a mass /r 2 ); must vanish in the test body limit 
when gravitational radiation vanishes, (i.e. be proportional to n); must be dissipative, 
or odd in velocities; must be linear in the spins; must be a vector, not a pseudovector; 
and must change sign under the interchange 1^2. Note that other possible terms, 
such as Lj\j( n ■ <S) can be seen to be linear combinations of the terms above using 
standard vector identities. The prefactor 1/5 is chosen for convenience. To make the 
terms of 0(e 7 / 2 ) beyond Newtonian order, Ag, Bg, Cg, A^, B^ and C| must be of O(e), 
and Dg and must be of 0(e 2 ). The only orbital variables available to construct 
expressions of the relevant order are v 2 , m/r and r 2 . Thus Ag, Bg, Cg and Dg can 
be written in terms of 15 arbitrary coefficients, in the form 

2 m .2 

,4,5 = a\v + a 2 h a 3 r , 

r 

2 m .2 

Bg = a^v + 05 h a^r , 

r 

2 m .2 

Cg — ajv + ag h agr , 

r 

2 

4 2-2 .4 2 m .2 m 171 <„\ 

Ms — a W v + anv r + ayir + a\^v h a^r h ais— o- . (9) 



In a parallel manner, we can write A^, B^, and in terms of its own set of 15 
coefficients. Because all expressions involving spin-orbit terms divide naturally into 
those involving the total spin S and those involving the spin parameter £, we can 
solve for each set using identical methods; we will focus on the S-terms. Our goal is to 
evaluate these thirty coefficients by imposing energy and angular momentum balance. 

Because the equations of motion at 2.5PN order and 3.5PN order have dissipative 
terms, the energy and angular momentum are no longer conserved explicitly. Further- 
more, they are ambiguous because one has the freedom to add arbitrary terms to E 
and J at 2.5PN order and 3.5PN order to redefine them without affecting their con- 
servation through 2PN order. Similarly, the spins are strictly defined only up to the 
order at which radiation reaction begins, and so one has the freedom to add a 3.5PN 
term to each spin, without changing its behavior at "conservative" orders. 

Adding to E and J the appropriate 2.5PN terms to account for the coefficients a 
and /3 in Eq. @, and adding the most general 3.5PN spin-orbit terms, with arbitrary 
coefficients, we can define new quantities E* , and J* according to 

E* — E N + SE 2 .5PN + 5-&J.5PN-SO , 

J* = M L N + 5 + £n x [n x (45 + 3£)] + <5J 2 . 5 PN + <5J 3 . 5 PN-SO - (10) 
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where, from our earlier work 12,13 , we can write 

SE 2 . 5FN = |a^r[(2 + a)« a - /?r 2 ] , 
5 r z 

xi n 8 n 2 m - 

dJ 2.5PN = q t — t^Ln, (11) 
5 r z 

and for the 3.5PN-SO expressions we write the general parametrized form 
<^3.5PN-SO = ~\^A * {(^N ' 5) ^Qit; 2 + a 2 r 2 + ^3^) + (S —> £) j , 
5J3.5PN-SO = "^^2 I 7 " 5 {j 1 "" 2 + f 27 " 2 +^3y) + rn(n ■ 5) ^ 

+v(n ■ 5) + 7 8 f 2 + 79-^-) + n ( v ■ S) (jiov 2 + 

+7l3rv(v-5) + (5-»€)l. (12) 



2 2 771 
1 74 V + 75r + 76 — 



.2 . m 
71 ir +712 — 
r 



where the notation S — * £ means repeat the preceding terms replacing S with £, with 
an appropriate set of arbitrary coefficients. This gives a total of 32 arbitrary coefficients. 

We now take time derivatives of E* and J* in Eqs. (|10p . substituting the Newtonian 
and PN spin-orbit accelerations explicitly, to obtain 

E* = flV ■ (a 2 .5PN + a 3 . 5 pN-S0) + -jT^-^.SPN + ^<^3.5PN-SO 1 

j* = <S3.5PN-SO + X (a 2 .5PN + a 3 . 5 pN-So) + ^^2.5PN + ^<5J3.5PN-So(13) 

In fact, we will show in Appendix [B] that, if we assume that (<Si)3.5PN_go is 
orthogonal to Si (and similarly for spin 2), then the most general 3.5PN expression 
for (5i)3.5pn-SO turns out to be a total time derivative, which can be absorbed into a 
meaningless 3.5PN correction to the definition of Si. Hence we can assume henceforth 

that S3.5PN-SO = (Sl)3.5PN-SO + (<S2)3.5PN-SO = 0. 

We now substitute the appropriate terms from the equations of motion ([3]) and {SJ , 
and calculate explicitly the time derivatives of the 2.5PN and 3.5PN-SO contributions 
to E* and J*. These time-derivative terms may be calculated using the identities 
shown in AppendixfS] which are derived using the Newtonian equations of motion and 
the 1PN spin-orbit terms. When evaluating d5E2,spa/dt and dS^2.5PN/di, in order 
to obtain all terms that contribute at 3.5PN-SO order, we must include the 1PN spin 
orbit terms present in the expressions in Appendix [X] The result is 



•» 8^i 2 m 2 2 n-2s 8^i 2 m f 

= ~ ) ~ Tip" Ln ' (4S + 3 ^ 



-(2 + a)v 2 +3f3f 2 



5r 5 \ 



(L N '5) [Viv +V2V r +V%r +V4V h V 5 r \-Vq— 5- 



(V-3r 2 + 2^) 



j* = _^L N (■>-<* - 2 < - m 
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&u ma I 1 ~ „. , , 

1 --^L N L N -(45 + 30+rnx 



5r 4 1 2r 2 

,2 



(v - -fn) x (45 + 3£) 



C" J c L 4.^ 2.2.^.4.^ 2m 2 m m 

br 6 V r r r z 



+n(n-5) VC7-U +7Csur + TZgr +TZiqv YH\\r Yk-\2—t 



+rn(v ■ 5) ^i 3 u 2 + TZ u r 2 + ^15™) + rv(n ■ 5) (fti 6 u 2 + ftiyr 2 + TZ W - 
+v(v ■ <S) (jl l9 v 2 + TZ 20 r 2 + n 2 i™) + (S -» £) j . (15) 

The first term in each of Eqs. (|14p and (|15p is the 2.5PN quadrupole, or Newtonian loss 
term, while the second term in each case comes from the spin-orbit correction terms in 
Appendix 1X1 applied to dSE2.5P^s/dt and gW^.spn/^- in the third set of terms in each 
case, the 27 coefficients V n , n — 1 ... 6 and TZ n , n = 1 ... 21 in the 5-dependent terms 
are functions of the 15 coefficients a n from the equations of motion ((8} and Q, and of 
the 16 coefficients a n and 7™ from the 3.5PN ambiguity terms in E* and J*. A parallel 
set of 27 coefficients appear in the ^-dependent terms, with identical dependences on 
the corresponding 15 + 16 arbitrary coefficients. 

We now use the assumption of energy and angular momentum balance to equate 
the rate of energy and angular momentum loss to the corresponding far-zone fluxes [2j 
[3|. The lowest-order Newtonian and the 1PN spin-orbit contributions are given by 



Efar 



zone 



2 _ 37t; 2 _ 12 m 



+(L N -£) (W 2 - I9v 2 -8™ 
• 8/i 2 m f 2 Q .2 . „ 

J tar zone = ~r o-L N (2u - 3r +2 

5 r A 



v 

m 



4 ^ 2 /^ ( R a- 



5 r 3 
5r 3 

+n(n-5) ^1 
+rn(v ■ <S) 



6n 



! ft 4 
— OU — 


50 2 m 

y w 7 + 


50 .2 m 


"2^ 


3 r 




— 30r 2 u 2 


+ 25, 2 ™- 
r 


„ .2 Wl 

f 6r — 

r 






- fv(n • 5) 


(l8i 2 - 


30r 2 4 



+v(v • <S) (6« 2 - 12r 2 + 23 — 



(16) 



in 

r 

^ e {-.A „ 2.2 10 4 22 2 m 23 . 2 m 

+5 5 r — iv r 1 1! 1 r =- 

V 3 3 r 3 r 3r z 

/ /.„ 4 „„.2 2 41 ?m „. 2 m 4m 2 \ 

+n(n ■ £) 13v - 20/ u H u h 6r 2 h — 5- 

Y 3 r r 2>r z J 

+rn(v ■ i) (lv 2 - 5r 2 - ^) - rv(n ■ £) (l3« 2 - 20f 2 + ^) 



s 



+v(v.e)(^ 2 -^ + ^)r 

After rewriting some of the terms in Eq. (115[l using standard vector identities, we 
compare Eqs. (|14p and (|15f) to Eqs. (|16|l and (|17p term by term to obtain 54 constraints 
on the 62 coefficients. It turns out, however, that 4 of these constraints are not linearly 
independent of others, so there are 50 non-trivial constraints, leaving 12 undetermined 
degrees of freedom. The specific choice of the free coefficients is somewhat arbitrary; 
one choice gives the following values for the coefficients ([9]) in the equations of motion 
©: 

a x = 2820[2160] + 1574 + 4577 + 4579 + 157n + 45712 - 3q 2 , 
a 2 = -1728[-1348] - I374 - 397 7 - 427 9 - II711 - 427 i2 + 3a 2 , 
+48[36](a-/3), 

a 3 = -6020[-4620] - 3574 - IO577 - IO579 - 357n - 1057i 2 + 7q 2 , 
a 4 = -220[-164] - 74 - 3 77 - 3 79 - 2 7u - 3 7l2 , 

a 5 = y [36] + 74 + 3 77 + 3 79 + 2711 + 3 7 i2 + 16[12]q , 
a 6 = 860[640] + 5 74 + 15 77 + 15 79 + IO711 + 15 7 i 2 , 
a 7 = -788[-608] - 4 74 - 677 - 157 9 - 4711 - 15712 , 
3152 

a 8 = -jj-[808] + 4 74 + I677 + 24 79 + 4 7 n + I6712 - 32[-24]a , 

a 9 = 2460(1900] + IO74 + 3O77 + 4579 + IO711 + 457i 2 , 

aio = -148 [-112] - 274 - 377 - 379 - 2711 - 3712 , 

an = 3320(2540] + 2574 + 6O77 + 6O79 + 25711 + 6O712 , 

ai2 = -6020[-4620] - 35 74 - IO577 - IO579 - 35m - 1057i 2 , 

1276(968] „.,-,„„ 
013 = h 474 + II77 + 979 + 4711 + 5712 , 

au = -4392[-3372] - 2374 - 8777 - 7879 - 23711 - 54712 + 48[36]a , 
872 " 

- — -274-877-679-2711-2712, (18) 



015 = —376 



where the numbers in square brackets represent the values to be used, along with the 
corresponding set of six free coefficients, for the terms in Eq. 1(8} involving £. 
The unique choice of the twelve coefficients 

45 287 89 140 263 
c*2 = — , 74 = — , 77 = — ;r, 79 = jp 711 = g", 712 = -1, (19) 

for the S terms, and 

105 181 155 Q , 105 fnn . 

Q-2 = g", 74 = — , 77 = g-, 79 = -34, 711 = — , 712 = 2, (20) 

for the { terms, along with the values a — —1, and (3 = for the harmonic Damour- 
Deruelle gauge at 2.5PN order, gives precisely the 3.5PN spin-orbit radiation reaction 
terms derived in 1101. 
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4 Gauge Freedom and Arbitrary Coefficients in the Equation of Motion 



The formulas for energy and angular momentum flux in the far zone are gauge invariant, 
while the equations of motion are gauge, or coordinate dependent. Any coordinate 
transformation x^ — > x^ + f, where is, in a suitable sense, of 2.5PN and 3.5PN 
order relative to x^, will induce changes in the variables of a binary system, such as 
the relative vector x and the spin vectors. Notice that a transformation of coordinate 
time simply induces a velocity-dependent change in x via x(i + St) = x(t) + vSt. 
As for the spin, any change induced by a gauge transformation at 2.5PN or 3.5PN 
order can always be reabsorbed into a new definition of spin, since it is ambiguous at 
radiation-reaction orders. Therefore we will only consider coordinate transformation 
induced changes in the relative vector x at 2.5PN and 3.5PN-SO orders, according to 



Sx 



2.5PN 



+ <5x 3 . 



5PN-SO ■ 



(21) 



The 2.5PN order coordinate change that corresponds to the arbitrary coefficients a 
and (3 in Eq. ((3} was calculated in [121113] . and is given by 



<5x 2 



5PN 



8fim 
15r 



[/3rn + (2(3 - 3a)v] 



We can derive directly 



(22) 



v = v + <5x 2 . 5 PN + <5x 3 . 5 p N _ S o , 
dv' _ dv 

~ dt 



df 
mx' 
r /3 



<5x 



2.5PN 



+ <5x 3 



5PN-SO • 



7T1X 171 

— + -3-(5x 2 . 5 PN - 3nn ■ <5x 2 . 5PN ) 
+-^(<5x 3 . 5 p N _ S o - 3nn ■ <5x 3 . 5PN -so) 



(23) 



The 2.5PN terms in these equations must also be used to determine the induced change 
in the 1PN spin-orbit acceleration terms in Eq. (5). In evaluating 5x 2 spn explicitly 
using Eq. (|22|l . the 1PN spin-orbit equations must be employed wherever an accelera- 
tion occurs. The result is that the equation of motion ^ changes between the original 
and the new coordinates by a quantity Q given by 



Q = 



8/ira 
5r 3 

8/im 
5r 5 



(zfiv 2 + (2a - 3/3)™ - 5/3f 2 ) rn - (v 2 - ™ - 3r 2 ) av 
N-SO - ^3 (5x 3 . 5 pn-so — 3nn ■ 5x 3 .5pn-so) 
^-L N • (45 + 3£) (3(a - (3)rn + av) - fv x (45 + 3£) (a + ~0 



in x (45 + 3£) (pv 2 - f3™ - (9a + 6/3)r 2 ) 



(24) 



Note that the 2.5PN terms in Eq. (|24p match exactly the arbitrary terms in Eq. (J3J). 
We now want to find a form for <5x 3 .5PN-SO so that the 3.5PN-SO terms in Eq. (12411 
match the terms in (|5J generated by the arbitrary coefficients in Eq. ©. This can 
be done either by direct integration to find 5x 3 .spn-SOi or by assuming a suitable 
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form for <5x3.5pn-SO an d seeing if one can solve for a set of coefficients. Remarkably, 
a solution can be found, and is given by 



<5x 3 .5P N _ so = -5^2 |~(^N ■ s ) (74 + 377 + 379 + 711 + 3712 - ~a 2 ) 
+ t^(Ln ■ S) ^74 + 377 + 379 + 3712 - ^ Q 2^ 

(77 + 79 + Ji2)v 2 + (74 + 377 + 379 + 711 + 37i 2 )r 2 



+ i nx5 



-(712-^)^ 
3 r 



-fv x 5(79 +712) + (S ->£) 



(25) 



The 12 (6 + 6) coefficients correspond precisely to the 12 degrees of freedom in Eqs. 
OH. 



5 Concluding remarks 



We have used energy and angular momentum balance to deduce the general form of 
the 3.5PN spin-orbit radiation reaction terms in the two-body equations of motion, 
and showed that the remaining undetermined degrees of freedom correspond to the 
freedom to change gauges or coordinates at the corresponding post-Newtonian order. 
A specific choice of the free coefficients yields 3.5PN spin-orbit terms in the equations 
of motion identical with those derived from first principles. The results were subject 
to the physically reasonable assumption that gravitational radiation reaction has no 
effect on the magnitude of the individual spins, to 3.5PN order. 

A natural extension of this work is to determine the contribution of spin-spin 
interactions in radiation reaction using balance arguments and to compare the results 
with those calculated from first principles by Wang and Will This work is in 

progress. 



A Extracting total time derivatives 

Using the Newtonian equations of motion plus the 1PN spin-orbit terms, it is straightforward 
to establish a number of identities, which may be used to relate collections of terms to total 
time derivatives of other expressions. For any non-negative integers s, p and q, we obtain 



d I v 2a rP 
dt \ ri 



d 
dt 



} 2s r P , 



v 2s-2j.p- 
r q + l 



4 1 1 \ 2 -2 n -2 2"" 
pv — (p + q)v r — 2sr pv — 



p v 
2 r 



— L N ■ (45 + 34) 



d I v 2a rP 



' dt 







r 3 1 




|.«( 


v rn 


{ r-9+2 




2 



(45 + 34)) 



(26) 
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Another set of identities, to be used only in 3.5PN terms, require only the Newtonian equations 
of motion: 



dt \ ri 



d ( v 2s fP 



dt \ ri 



d ( v 2s rP 



dt \ ri 



,,2s-2^,p-l 



,.q- 



.p-i (■ 
1 I 



pv — (p + q)v r — 2sr 



+2v 2 rrx < - i v j) } , 
,,2s-2 f p-l 



-2ra^-x^v^ 



J ~ ^+1 \ 



pv — (p + q)v r — 2sr 



■ pv 



pv — (p + q)v r — 2sr pv — 



+v 2 rr ^v l v-i n l n J ^ ^ 



- pv 



(27) 



B Evolution of spins at 3.5PN order 



In this appendix we justify our assumption that the individual spins are unaffected by 3.5PN 
spin-orbit effects, i.e. that (dSi/dt)^ 5 pj^_ go = 0, and similarly for body 2. First we write 
down the general form that 3.5PN spin-orbit terms could take, consistent with the assumptions 
used in earlier sections, namely 



A/iSi + M>nn • Si + Af 3 rnv ■ Si + A/" 4 r-vn ■ Si + A/5 vv • Si > , (28) 



(Sl)3.5PN-SO 



where Ai and A2 are each linear combinations of v* , u 2 r 2 and v 2 m/r etc. at 0(e 2 ) (containing 
6 terms each), and A/3, A/4 and A/5 are each linear combinations of v 2 , r 2 and m/r. Note that 
other possible terms, such as Ln X 5i, or LjmLn • Si can be rewritten as linear combinations 
of the terms above. 

We now impose the physically reasonable constraint 1S1 ■ (5i)3.spn— so = 0, which implies 
that radiation reaction does not change the magnitude of the body's spin, only its orientation. 
That constraint implies that only the third and fourth terms in Eq. I|28j l survive, and then 
only in an antisymmetric combination that leaves (<Si)3.5pn— SO m the general form 



(5i) 3 



-rLjv x Si 



(civ 2 



c 2 r 



c 3 - 



(29) 



However, it is straightforward to show, using the identities in AppendixfS] that the right-hand 
side of Eq. H29H can be written as a total time derivative and therefore can be absorbed into 
Si, independently of the values of ci, C2 and C3. This is in accord with the result derived from 
first principles in 1101 . 
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